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ABSTRACT 

We report on the observation of the region around supernova remnant G65. 1+0.6 with the 
stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL 
J1954.3+2836 and 1FGL J1958. 6+2845. They are identified as GeV pulsars and both have a 
possible counterpart detected at about 35TeV by the Milagro observatory. MAGIC collected 
25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. 
We therefore report differential flux upper limits, assuming the emission to be point-like (< 0.1°) 
or within a radius of 0.3°. In the point-like scenario, the flux limits around 1 TeV are at the level 
of 3% and 2% of the Crab Nebula flux, for the two sources respectively. This implies that the 
Milagro emission is either extended over a much larger area than our point spread function, or 
it must be peaked at energies beyond ITeV, resulting in a photon index harder than 2.2 in the 
TeV band. 

Subject headings: pulsars: individual(PSR J1957+2831, LAT PSR J1954+2836, LAT PSR J1958+2846) 
— ISM: supernova remnants 1 



1. Introduction 

In February 2009, the Fermi collaboration pub- 
lished a list of the most significant gamma-ray 
sources above 100 MeV, detected by the large 
area tel escope (LAT) with in 3 months of obser- 
vation (jAbdo et al. 2009b ). The energy spec- 
tra of these sources often extend to several GeV, 
where at some point the steeply falling flux lev- 
els are too low to be detected by the limited de- 
tection area of a satellite instrument. Many of 
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the LAT sources are hosted by our galaxy, and 
34 of those are within the field of view of the 
Milagro gamma-ray observatory, wh ich was lo- 



cated n ear Los Alamos, New Mexico (jAbdo et al 
(|2009al ) and references within). The Milagro col- 
laboration, therefore, reinvestigated their previ- 
ously gained skym ap to look for cou nterparts to 
these GeV sources ( Abdo et al.ll2009a ). The sensi- 
tivity of the instrument peaks in the energy range 
10 — 50TeV, although it ultimately depends on 
the energy spectrum and the declination of the 
source. Due to the reduced trial factor, 14 new 
Milagro sources could be claimed with confidence 
levels above 3 a. 

Two of these Milagro-detected Fermi bright 
sources are in the vicinity of G65. 1+0.6, a 
faint supernova r emnan t ( SNR) first repo r ted by 
lLandecker et~ai] (|l990l ). iTian fc Leahvl (|2006h 



suggested its distance to be 9.2 kpc and a Se- 
dov age of 40 — 140 kyr. An association to the 
radio pulsar PSR J 1957+2831 was suggested by 



Lorimer et all (|l998h . 



The gamma-ray emission in the region of 
G65 1+0.6 was first detected by the COS-B satel- 
lite (jSwanenburg et al.lll98ll ) as 2CG065+00, and 
later confirmed by the EGRET sate llite (3EG 
J1958+2909) in Irlartman et~ai1 (|l999h . where a 
possible extension or multiple sources were de- 
noted. As of now, the two sources could be de- 
tected as two individual sources by the LAT as 
1FGL J1954.3+2836 and 1FGL J1958.6+2845, as 
reported in the first year catalog of Fermi sources 
( Abdo et al. 2010bl ). They were analysed and re- 
ported as gamma-ray pulsars found through blind 
search in Saz Parkinson et al. (|2010() . Abdo et al 
(|2009ch and lAbdo et alj (|2010al ) Their peri- 
ods (290ms, 92.7ms), spin-down luminosities 



(104.8 x 10 34 ergss-\ 33.9 x lO^ergss" 1 ) and 
energy cutoffs (2.9 GeV, 1.2 GeV), but also the 
characteristic magnetic fields and ages (69.5 kyrs, 
21 kyrs) lie in the average range of all Fermi pul- 
sars. In the following, we will use the names J1954 
and J1958 for these 1-year catalog sources, and 
J1954n/J1958n to specifically refer to the 3-month 
bright source list positions. 

The detections by Milagro revealed signifi- 
cances of 4.3 a for J1954 and 4.0 a for J1958 
(|Abdo et al.ll2009ah . Flux values are stated for a 
characteristic median energy of 35 TeV. The an- 
gular resolution of Milagro is about 0.4 — 1.0°, so 
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these values can be expected to hold also for the 
1-year catalog positions of the sources, which are 
offset by < 0.1°. Since gamma-ray pulsars typi- 
cally have energy cutoffs as low as several GeV, 
the Milagro signals, if real, can be expected not 
to be caused directly by the pulsars, but possi- 
bly by associated objects, such as a pulsar wind 
nebula (PWN) or, in the case of J1954, the shell 
of the SNR, which surrounds the pulsar. Gamma 
rays at TeV energies may also be produced in 
an interaction of the shell with a possibly coin- 
cident molecular cloud. This cloud may be lo- 
cated at the position of the infra-red source IRAS 
19520+2759, which has associated CO line, H 2 
and OH maser emission at a similar distance as 



the SNR (lArauilfa fc Kwoklll987ft . 

The MAGIC telescopes use the Cherenkov 
imaging technique and are located on the Ca- 
nary Island of La Palma (28.8°N, 17.8°W, 2220 
m a.s.L). It is the instrument with the lowest 
energy threshold among all Cherenkov telescopes. 
In single-telescope observations, as presented here, 
the nominal thre shold in low zenith a ngle observa- 
tions is 60 GeV (| Albert et al.ll2008ah . It is there- 
fore the instrument of choice to connect the upper 
ends of the Fermi spectra, which typically end at 
some tens of GeV, with the detections provided 
by Milagro at some tens of TeV. To bridge these 
spectra in the TeV range, and identify possible 
object associations for the Milagro signals were 
the main motivations for our investigation. 

2. Observations 

Motivated by the fact that J1954o was not 
marked as a pulsar in the initial Fermi Bright 
Source List, we observed J1954o as a main target 
in July and August 2009. The observations were 
carried out in false source tracking (wobble) mode 



(|Fomin et al.l 119941 ) . which yielded two datasets 
with offsets of ±0.4° in RA from this source, see 
Figure [TJ The wobble position was altered every 
20 min, and the data were taken at zenith angles 
between and 43 degrees. At the time, the sec- 
ond MAGIC telescope was still under comission- 
ing, so the analysis presented here used only the 
data from the stand-alone MAGIC-I telescope. 

Quality selection cuts were applied to the event 
rate, the spread of hardware-sensitive shower im- 
age parameters, and few parameters that charac- 



terize the transparency of the atmosphere, such as 
the sky temperature and humidity. After all data 
selection cuts, 25.5 hours of high quality data were 
left for the analysis of J1954. 

Besides this, we took advantage of the fact that 
J1958 is in the field of view of one of the two wob- 
ble positions. It could be analysed in a specific 
way described in the next section, yielding 13.8 
hours of effective observation time. 

3. Analysis and Results 

The data w ere analysed in the MARS analy- 
sis framework ( Moraleio et al.l l2009). which is the 
standard software used in the analysis of MAGIC 
data. After the air shower images of the photomul- 
tiplier tube camera are calibrated, and times and 
charges of each pixel are extracted, a three-stage 
image cleaning is applied to filter out uncorrelated 
noise from the data acq uisition electroni cs and the 
night sky background ( Aliu et al. 2009). Shower 
image paramet ers are calculated, an d the Random 
Forest method ( Albert et al.ll2008bT) is used to de- 
rive estimators for the shower direction, its energy, 
and its likeliness to be of hadronic origin. 

The observational setup described in the previ- 
ous section requires a different analysis treatment 
of both sources. The main target of the observa- 
tion, J1954, can be analysed using standard wob- 
ble analysis procedures. This means that the pho- 
ton flux from the source is compared to the one on 
the opposite side of the camera (anii-source, solid 
squares in Figure [T]) . In this way, exposure inho- 
mogeneities that can arise from imperfections in 
the photo multiplier tubes, trigger electronics or 
the signal transmission, cancel out, because both 
wobble positions were equally populated. 

In the case of J1958, a wobble analysis, us- 
ing only one of the two wobble positions, is not 
guaranteed to cancel out these inhomogeneities. 
Therefore, the analysis was done in ON/OFF 
manner, using the near wobble sample as ON- 
source data and the far wobble sample as OFF- 
source data (hollow squares in Figure [l}. Having 
the OFF-source at the same position in relative 
camera coordinates as the source in the ON sam- 
ple, the exposure inhomogeneities cancel out. 

To test the presence of a gamma-ray signal, the 
distributions of squared angular distances (9 2 ) be- 
tween photon directions and the source positions 
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were used. In these 6* 2 -plots, a signal is expected to 
produce an excess where 9 2 approaches zero. How- 
ever, the integrals over the expected signal regions 
of the 9 2 distributions agree, for both sources, well 
with the corresponding integrals done with respect 
to the OFF regions. Furthermore, the shapes of 
the ON and OFF distributions agree well with 
each other, and are sufficiently flat to exclude the 
unlikely case of an emission occuring by chance at 
both ON and OFF locations at a similar flux level, 
supressing the significance of the 9 2 comparison. 
Such a coincidence was also excluded by cross- 
checking the analysis using several OFF regions, 
and by thoroughly investigating the skymaps with 
different background estimation algorithms. Con- 
sequently, we report the absence of a significant 
signal for both sources. 

Besides this analysis, which takes advantage of 
a-priori defined source locations, also a 2.8° x 2.0° 
skymap of the area was investigated at different 
energy ranges. The trial factors implied by these 
searches were typically between 120 and 260, de- 
pending on the PSF, which is smaller for high en- 
ergies. Taking into account these trials, no signif- 
icant signal was found at any energy. 

We convert our data into three differential flux 
upper limits for each source. To do so, the 
data was divided into three bins of estimated en- 
ergy, delimited by 120 GeV, 375 GeV, 2.8 TeV and 
12 TeV. For each bin, an event number upper limit 
is calculated from the above-mentio ned fl 2 -plots 
at 95 % confidence level (c.l.) after iRolke et al 



(|2005l ). and assuming an efficie ncy systematic er- 
ror of 30% (| Albert et al.ll2008al ). For each limit, a 
power law energy spectrum with a photon index 
of —2 is assumed both for the simulation of the 
effective area, and the conversion of event number 
upper limits to flux upper limits. The influence of 
this photon index is minor, since the energy ranges 
are sufficiently small. 

Finally, since the data can only be selected by 
estimated energy, a Monte Carlo simulation (MC) 
is used to estimate the median true energy of the 
data that remains after all cuts. 

This analysis assumes a source extension sim- 
ilar or smaller than the point spread function 
(PSF) of MAGIC. In an identically conducted 
analysis of contemporary Crab Nebula data the 
width of this PSF (defined as the Sigma of a two 
dimensional Gaussian function) was found to be 



about 0.08°. Since the Milagro source may be 
extended, a second analysis was done, assuming 
an extension of 0.3° instead. Since this is done 
simply by increasing the signal integration radius, 
more background events are included, which leads 
to higher upper limits. 

4. Discussion 

The derived 95% c.l. flux upper limits are sum- 
marized in Table [TJ Figures [5] and [3] display them 
in the context of the published satellite data and 
the Milagro flux estimation. The differences be- 
tween the limits of J1954 and J1958 are all com- 
patible with the statistical fluctuations that can 
be expected for 95% upper limits. Around 1 TeV, 
where MAGIC is most sensitive, the flux is lim- 
ited to 3% of the Crab Nebula flux for J1954 and 
2% for J1958. Assuming that the Milagro emis- 
sions originate from objects spatially coinciding 
with the Fermi pulsars within our PSF, the pho- 
ton index in the energy range of 1 to 35 TeV must 
be harder than 2.2 for J1954, and 2.1 for J1958. 
The spectral energy distribution is thus likely to 
peak at energies in excess of 1 TeV. 

If an extension up to 0.3° is assumed, the corre- 
sponding flux limits in Crab Nebula units are 14% 
for J 1954 and 3% for J 1958. In this extended case, 
the photon indices are limited to < 2.6 and < 2.2, 
respectively. It shall be noted that the biggest 
TeV pulsar wind nebulae have sizes of few tens of 
parsecs, which at the distance of G65. 1+0.6 would 
be within these 0.3°. 

Bringing together the existing flux data and 
our upper limits, we conclude that the most 
likely scenario to explain the gamma ray pro- 
duction measured by Milagro might be the ex- 
istence of two PWN, associated with the Fermi 
pulsars. With the ages of the pulsars being 
69.5 kyrs and 21kyrs, it is reasonable to ex- 
pect an inverse compton component that dom- 
inates their ene rgy outflows and may addition- 
ally be extended (Ide Jager &: Diannati-Atai'T 2009; 



Tanaka &: Takahara 2009f) . Such old, extended 
PWN are common TeV gamma-ray sources and 
frequently have an emission spectrum that peaks 
at TeV energies or above. 
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5. Summary 



We took 25.5 hours of good quality data in the 
area of the faint supernova remnant G65.1+0.6. 
In that region, the Milagro collaboration reported 
the emission of gamma rays with median en- 
ergy of 35 TeV in the vicinity of the two GeV 
Fermi pulsars 1FGL J1954.3+2836 and 1FGL 
J1958. 6+2845. Our observations, which were 
aimed to locate the source of the Milagro emis- 
sion, yielded no significant gamma-ray signal for 
these two a-priori known source locations. Also, 
no post-trial significant signal from a skymap of 
the area could be established. We extracted three 
differential flux upper limits for each source, as- 
suming two different extension radii. They arc 
summarized in Figures [2] and [3] 

Assuming that the two 4 and 4.3<r detections 
of Milagro are not statistical fluctuations, but real 
signals, the flux upper limits support the scenario 
in which the multi-TeV emission measured by Mi- 
lagro is caused by a different mechanism or object 
than the Fermi emission. Given the ages of the 
pulsars and the SNR, the existence of two very 
old pulsar wind nebulae, powered by the two GeV 
pulsars, seems very likely. 
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Gemeinschaft, and by grant DO02-353 of the the 
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Fig. 1. — Observation setup for the two Fermi 
sources J1954 and J1958 in the context of SNR 
G65.1+0.6 and a Milagro significance contour. 
J1958 appears only in one wobble position (Wl), 
so the OFF data is taken from the other wobble 
sample, using the same position relative to the 
pointing direction. The outline of the remnant 



i s tak en from the radio map in lLandecker et al 
(1990). The extension of the Milagro significance 



contour ( Abdo et al. 2009al ) is compatible with 
their point spread function. 
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Fig. 3. — Compilation of flux measure- 
ments and upper limits for 1FG L J195 8.6+2845 
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Fig. 2. — Compilation of flux measurements and 
upper limits for 1FGL J1954.3+2836 from Ferm i 



( Saz Parkinson et al. 2010 : Abdo et all 



2010M ). 



MAGIC and Milagro (|Abdo et al.M2009al) . The 
3% fraction of th e MAGIC Crab spectrum 



([Albert et al.ll2008al ) is shown for comparison. 
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Table 1 
Differential Upper limits 



Source Name 


Assumed Extension 


^mcd 


Significance 


F 95% 
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(GoV) 




(10~ 12 TeV _1 cm- 


- 2 .- 1 ) 


(10~ 12 TcVcm- 2 S - 1 ) 
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< 0.08° 
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-1.8 




27 
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+ 1.9 
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